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Abstract 



Jet fragmentation in pp and PbPb collisions at a centre-of-mass energy of 2.76 TeV per 
nucleon pair was studied using data collected with the CMS detector at the LHC. 
Fragmentation functions are constructed using charged-particle tracks with trans- 
verse momenta px > 4 GeV/c for dijet events with a leading jet of Pt > lOOGeV/c. The 
fragmentation functions in PbPb events are compared to those in pp data as a function 
of collision centrality, as well as dijet-px imbalance. Special emphasis is placed on the 
most central PbPb events including dijets with unbalanced momentum, indicative of 
energy loss of the hard scattered parent partons. The fragmentation patterns for both 
the leading and subleading jets in PbPb collisions agree with those seen in pp data 
at 2.76 TeV. The results provide evidence that, despite the large parton energy loss 
observed in PbPb collisions, the partition of the remaining momentum within the jet 
cone into high-px particles is not strongly modified in comparison to that observed 
for jets in vacuum. 
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1 Introduction 

When colliding lead nuclei (PbPb) at a nucleon-nucleon centre-of-mass energy of ^/s^ = 
2.76 TeV, one expects to form a system of hot and dense matter at energy densities that have 
not been explored before. One of the early proposed experimental signatures of the formation 
of a dense system in such collisions was "jet quenching", i.e. the suppression or disappearance 
of the spray of hadrons resulting from the fragmentation of a hard scattered parton having 
suffered energy loss in the medium |IT|. The energy lost by a parton in the produced medium 
provides fundamental information on its thermodynamical and transport properties |21|3|. Re- 
sults on the suppression of inclusive hadron production at high transverse momenta (pj), as 
well as on the modified high-px dihadron angular correlations obtained from nucleus-nucleus 
collisions at ^/s^ = 200 GeV at the Relativistic Heavy Ion Collider (RHIC) [4 - 7J have shown 
the existence of partonic energy loss in dense strongly interacting matter. Similar observations 
have been made at the Large Hadron Collider (LHC) [.81^14]. 

At LHC energies, high-px jets have been fully reconstructed in heavy-ion collisions. A sig- 
nificant dijet transverse momentum imbalance is observed, when comparing to a reference 
distribution corresponding to pp collisions at the same nucleon-nucleon centre-of-mass en- 
ergy llT5l - [T7| . Such an observation is consistent with energy loss of the hard scattered partons 
in the dense medium produced in central PbPb collisions. In the same set of results, the re- 
distribution of the lost jet energy is studied using jet-track correlations [16J. It is found that 
the missing pj opposite to the leading jet can be recovered only by summing up the contribu- 
tions of particles down to a low pj of 500 MeV/c with respect to the beam axis and out to large 
radii in pseudorapidity and azimuthal angle, AR = (Ac/))^ + [Ar])^ > 0.8, with respect to 
the jet axis IIT6l . Since these results show that the quenched energy is transferred out of the jet 
cone, the jet clustering algorithm reconstructs jets with a reduced energy from the fragments 
of the partons after they have lost energy in the medium. The study presented here inves- 
tigates to what extent the fragmentation pattern of partons that have traversed the medium 
resembles vacuum fragmentation, by constructing fragmentation functions in PbPb collisions 
and comparing them to those from unquenched jets, as measured in pp collisions. Fragmen- 
tation functions encode the probability for a parton to fragment into particles carrying a given 
fraction of the parton energy. Colour-charged partons undergo showering processes into par- 
tons of successively lower energy which hadronise into colour-neutral final-state particles. The 
evolution of such a parton radiation and splitting process leads to a characteristic shape of the 
fragmentation function [18|. Theoretical models of jet quenching predict an effective change of 
the shape of the fragmentation function due to the change of the parton radiation pattern in the 
medium [19-22J. Experimentally, fragmentation functions are constructed by correlating the 
reconstructed jet momentum with the momenta of charged particles projected onto the jet axis. 
The jets are defined using the final-state particles produced in the collision, clustered with the 
anti-A;x jet algorithm [23J. In this Letter, we present a measurement of fragmentation functions 
in pp and PbPb collisions and a detailed comparison of their shapes measured in the two sys- 
tems. The measurement is restricted to the high-px component of the fragmentation function, 
using charged particles of px > 4GeV/c that lie within AR < 0.3 around the reconstructed jets. 

2 Experimental Setup 

The Compact Muon Solenoid (CMS) detector is described in Ref. p4) . Only the detector sys- 
tems used in this analysis are discussed hereafter. The central part of the CMS detector contains 
a superconducting solenoid that provides a homogeneous magnetic field of 3.8 T parallel to the 
beam axis. Charged-particle trajectories are measured using silicon pixel and strip trackers that 
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4 Jet and Track Reconstruction 



cover the pseudorapidity region \ri\ < 2.5, where t] is defined as rj = — log [tan {9/2)] and 9 
is the polar angle with respect to the anticlockwise beam direction. An electromagnetic crystal 
calorimeter (ECAL) and a brass/ scintillator hadron calorimeter (HCAL) surround the tracking 
volume and cover < 3.0. The ECAL calorimeter is segmented in quasi-projective cells of a 
granularity in pseudorapidity and azimuthal angle of At] x = 0.0174 x 0.0174 in the barrel 
(I //I < 1.5), increasing across the endcap (1.5 < \rj\ < 3.0) to 0.09 x 0.09 at \f]\ = 3.0. The 
HCAL has a cell granularity of A?/ x Ai^ = 0.087 x 0.087 in the barrel region and a variable cell 
granularity, changing as a function of in the endcap region [24 J. A forward steel/ quartz-fibre 
Cherenkov hadron calorimeter (HE) extends the coverage to |f/| = 5.2. The CMS trigger sys- 
tem is composed of a first level made of custom hardware processors, which use information 
from the calorimeters and muon detectors to select events, and the High-Level Trigger (HLT) 
processor farm, that further decreases the event rate, before data storage. 

3 Data Selection 

The PbPb and pp data analysed in this Letter were collected with the CMS detector in 2010 
and 2011, respectively, at a centre-of-mass energy of 2.76 TeV per nucleon pair. The integrated 
luminosities for the PbPb and pp data samples used for this analysis are Lint ~ 6.8 fA)^^ and 
Lint ~ 231 nb^^, respectively. The HLT system is used to select events containing high-px 
jets reconstructed from calorimeter towers. In PbPb collisions, the trigger threshold is pj = 
35GeV/c, applied on the raw calorimetric jet energy. Eor pp collisions, events are selected if 
they pass a jet trigger threshold of pj = 40GeV/c on the calorimetric jet energy. As found 
in Ref. 1161 , for the jet selection used in this analysis, requiring a 100 GeV/c jet in |?/| < 2, 
the triggers are more than 99% efficient. In addition to the trigger decision, standard event 
selection criteria are applied I16J , including a rejection of beam related backgrounds, a selection 
of inelastic hadronic collisions by requiring a two-sided coincidence of signals in the HE and a 
well-reconstructed event vertex. 

Eor the analysis of PbPb data, it is important to determine the degree of overlap between the 
two colliding nuclei in each event, termed collision centrality. Centrality is determined using 
the sum of transverse energy reconstructed in the HE. The distribution of the HE energy is used 
to divide the event sample into percentiles of the total nucleus-nucleus interaction cross section. 
Eor the purpose of this analysis, the data set is split into two centrality bins, the 0-30% most 
central events (i.e. those which have the largest overlap between the two colliding Pb nuclei) 
and the remaining peripheral events in the 30-100% centrality range. A detailed description of 
the centrality determination can be found in ||T6ll . 

4 Jet and Track Reconstruction 

Eor both pp and PbPb collisions, the analysis is based on jets reconstructed using the anti-Zcj jet 
algorithm, with a radius parameter (R) of 0.3, utilizing particle-flow (PE) objects that combine 
tracking and calorimetric information Il25ll26ll . In the PbPb data, the contribution of the under- 
lying heavy-ion event is removed using an iterative pileup subtraction method [27J . Since this 
procedure determines the underlying-event background using data outside the jet, the result is 
insensitive to the fragmentation properties of the jet. The jet-finding efficiency is above 95% for 
jets of Pt > 40 GeV/c, and above 99% for jets of px > 50 GeV/c ||25| . The relative jet momentum 
resolution in pp collisions is found to be 19 (13)% at px = 40 (100) GeV/c, improving with jet 
momentum. In central PbPb collisions, the momentum resolution deteriorates to 24 (16)% at 
Px = 40 (100) GeV/c II25II due to fluctuations in the underlying event. Eor both pp and PbPb 
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events, the jet momentum response has little or no deviations from a Gaussian shape. 

In all cases the reconstructed jet momenta are corrected to final-state stable particle (lifetime 
T with CT > 10 mm) level using factors derived from PYTHIA 6.422 l28] (tune D6T [291 130| , 
CTEQ6L1 PDFs [31J) pp simulations at = 2.76 TeV [32J. The uncertainty in the corrected 
jet energy scale is about 3% for pp events, resulting in a per-bin jet-yield uncertainty of ±15%. 
In the case of PbPb events, due to the influence of the underlying event, the uncertainty in the 
jet energy scale increases to about 4% for peripheral events (30-100% centrality) and 5% for 
central events (0-30% centrality) which results in per-bin jet-yield uncertainties of ±20% and 
±25%, respectively. 

The dijets selected for this analysis consist of a leading jet (denoted by subscript 1) with px,i > 
100 GeV/ c and a subleading jet (subscript 2) of Pt,2 > 40 GeV/ c, with axes that lie within \rj\ < 2. 
The Pi thresholds are chosen to ensure high reconstruction efficiency for the leading and, espe- 
cially, the subleading jet. In addition, the azimuthal opening angle A(^i 2 between the leading 
and subleading jet is required to be larger than In/ 3. No explicit requirement is made on the 
presence or absence of a third jet in the event. 

A detailed description of the charged-particle reconstruction algorithm and its performance 
can be found in Ref. El. The track-finding efficiency in the kinematic range of this study is (60- 
70)% and the corresponding correction is applied as a function of track pj, jet pi, and event 
centrality by reweighting the found tracks with the inverse of the reconstruction efficiency. 
The track reconstruction efficiency correction is derived from a Geant4 [33 | simulation of the 
CMS detector applied to pythia events, which are embedded into PbPb collisions simulated 
using HYDJET [34 1 in order to include the effect of the underlying PbPb event. The momentum 
resolution of the track reconstruction is a{pj)/ pj ^ 1-3%. 

CMS, pp, Vs = 2.76 TeV, L^^ = 231 nb"' CMS, PbPb, Vs^ = 2.76 TeV, L^^ = 6.8 nb"' 
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Figure 1: Dijet asymmetry, Aj, distributions in (a) pp collisions, (b) peripheral (30-100%) PbPb, 
and (c) central (0-30%) PbPb collisions. Data are shown as black points while the histograms 
show PYTHIA dijets, which when compared to PbPb data have been embedded into hydjet 
events. Error bars represent the statistical uncertainty. 

The dijet momentum balance is studied in terms of the dijet asymmetry ratio |[T5tiT7| , 

Aj = FiA^m^ (1) 

Pt,i ± Pt,2 

which is positive by construction. Figure [l] shows the Aj distributions in (a) pp and in (b,c) 
PbPb for two bins in event centrality. Central PbPb events (0-30%) show a significant excess of 
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unbalanced pairs when compared to both peripheral PbPb collisions (30-100%) and pp data. 
This can be interpreted as a direct observation of parton energy loss in central PbPb collisions. 

5 Fragmentation Functions 

The fragmentation functions are measured by correlating reconstructed charged-particle tracks 
falling within the jet cones, with the axis of the respective jet l|35ll . As done in previous mea- 
surements at hadron colliders [36 1, the fragmentation function is presented as a function of the 
variable 

p track 

where pj^^'^'^ is the momentum component of the track along the jet axis, and pjet is the mag- 
nitude of the jet momentum within the jet cone. The momentum components and the angle 
between the charged-particle and the jet axis are calculated in the dijet centre-of-mass frame, 
obtained by an approximate Lorentz transformation along the beam axis in the form of a pseu- 
dorapidity shift, defined as ?/dijet = (^/i + '/2)/2. The tracks are selected to lie within a cone of 
AR < 0.3 around the jet axis. The fragmentation functions, defined as l/Njet dN^^ack/ <^^, are 
normalised to the total number of selected leading or sub leading jets, Njet, respectively. To min- 
imise the contribution of tracks from the underlying event, only tracks with pj^'^^ > 4 GeV/c are 
selected. This restricts the measurement of the fragmentation function to the region of ^ < 4.5. 
The remaining underlying event contribution, not associated with the jet, is estimated by se- 
lecting tracks that lie in a background cone, obtained by reflecting the original jet cone about rj 
= 0, while keeping the same (p coordinate. The background contribution is accumulated jet-by- 
jet over the full event sample and subtracted to obtain the final dN^^^^^f distribution. Due 
to this procedure, jets in the region \rj\ < 0.3 are excluded from the analysis, to avoid overlap 
between the signal jet region and the region used for background estimation. 

Figure |2] shows the reconstructed leading and subleading jet fragmentation functions in pp 
collisions (right panel) and the corresponding jet px distributions (left panel), to illustrate the 
kinematic range in which the fragmentation functions are measured. Note that the higher jet 
momentum of the leading jet compared to the subleading jet, leads to an increased number of 
particles passing the pj'^'^^ > 4 GeV/ c selection for the fragmentation function measurement. 
This results in the observed excess of dNtrack/ d^ at high values of ^ for the leading jet over the 
corresponding distribution for the subleading jet. The pj^'^^ threshold on the tracks introduces 
a jet-px-dependent kinematic limit in the high ^ part of the spectrum. For a direct comparison 
between pp and PbPb collisions, the jet momentum resolution deterioration in PbPb events has 
to be taken into account. For this purpose, the reconstructed px of every jet in the pp data is 
smeared by the quadratic difference of the underlying event fluctuations in PbPb and pp. Fur- 
thermore, in order to keep the kinematic constraints consistent, a jet-px-dependent reweighting 
is applied to the pp data, after fluctuation smearing, so that the resulting jet px distribution 
matches that in PbPb. The reweighting factor is applied to each jet when generating the frag- 
mentation function for pp. The pp-based reference distributions obtained this way ensure that 
the jet fragmentation functions in PbPb and pp are compared for matching jet pj spectra. 

Figure |3] shows the fragmentation functions for (a) peripheral and (b) central PbPb collisions, 
for both the leading and subleading jets, compared to the pp-based reference. The ratios be- 
tween the PbPb fragmentation functions and the pp-based reference distributions are shown 
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CMS, pp, is = 2.76 TeV, L, = 231 nb ' 



CMS, pp, lis = 2.76 TeV, L = 231 nb"' 
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Figure 2: Data from pp collisions. Left: Leading and subleading jet pj distributions (not cor- 
rected for jet-finding efficiency and not unfolded for the jet energy resolution). Right: Frag- 
mentation functions reconstructed for the leading (open circles) and subleading (solid points) 
jets. The statistical uncertainties, shown as error bars, are smaller than the symbols in most 
cases. 

in panels (c) and (d). The corresponding jet px distributions illustrating the kinematic range 
of the measurement are shown in panels (e) and (f). The overlaid histograms in the same set 
of figures show the pp-based reference distributions. The systematic uncertainty, represented 
by the boxes at each point in panels (c) and (d), is obtained from the propagated jet and track 
reconstruction imcertainties. These comparison plots show that the shape of the fragmentation 
functions in pp and PbPb collisions agrees within uncertainties at all centralities for the leading, 
as well as for the subleading, jets. 

The uncertainties in the jet response may affect the results in different ways: smearing of jet 
energy due to fluctuations distorts the observed fragmentation functions, a miscalibration of 
the overall energy scale shifts the fragmentation function along the ^ axis, and a residual offset 
in the jet energy introduces a tilt of the shape of the distribution. These effects are studied using 
a Monte Carlo (MC) simulation, by varying the corresponding generator-level jet properties 
within the limits of the jet response uncertainty. The systematic imcertainties are determined 
by comparing the resulting fragmentation functions in the modified sample to the original MC 
reference. 

The systematic imcertainty due to the charged-particle reconstruction efficiency is obtained by 
comparing fragmentation functions based on efficiency-corrected tracks, with the fragmenta- 
tion functions using the MC generator information. Since the particle-flow event reconstruction 
algorithm uses reconstructed charged-particles for the jet pi determination, a failure to recon- 
struct a high-px charged-particle can lead to an underestimation of the jet momentum, resulting 
in an artificially high Aj measurement. The modification of the fragmentation function mea- 
surement due to this effect is studied in pythia + hydjet events. Since the Aj distributions in 
data and simulation are different, the magnitude of the corresponding effect in the PbPb data 
is estimated based on the reconstructed Aj distributions, and is accounted for in the combined 
systematic uncertainty. The effect of momentum resolution on reconstructed charged-particle 
tracks is estimated by smearing pythia + hydjet generator-level information. This is found to 
have little effect on the fragmentation function, in comparison to the unsmeared generator level 
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CMS, PbPb, ^ = 2.76 TeV. L = 6.8 |.ib^ 
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Figure 3: (a,b) Fragmentation functions reconstructed in peripheral and central PbPb data for 
the leading (open circles) and subleading (solid points) jets. (c,d) Ratio of each PbPb fragmen- 
tation function to its pp-based reference. Error bars are statistical, the hollow boxes represent 
the systematic uncertainty for the leading jet, and gray boxes show the systematic uncertainty 
for the subleading jet. (e,f) Jet px distributions in PbPb data (not corrected for efficiency and 
not unfolded for pj resolution) compared to the pp-based reference (see text). Only statistical 
uncertainties are shown in panels a, b, e and f . 



7 



information. The above uncertainties are combined in quadrature to give the total systematic 
uncertainty. 

Another potential source of systematic uncertainty comes from the response of the PF technique 
to jets with very different fragmentation functions. Since one component of the jet energy cor- 
rection accounts for the loss of low-px particles, a large change in the contribution of such tracks 
will result in an incorrect reconstructed jet energy. To study this effect, we have compared the 
response to separate quark and gluon jets from PYTHIA, whose fragmentation functions differ 
by 20-40% in the region 2 < ^ < 4 and by larger factors for ^ < 1. These dramatic differ- 
ences in fragmentation pattern do result in systematic differences in the jet momentum but the 
effect is only at the few percent level l|25| . The change in the reconstructed fragmentation func- 
tion resulting from the tiny shifts in the ^ parameter caused by these jet momentum offsets are 
negligible. 

To study in more detail the potential effect of medium-induced energy loss on the fragmenta- 
tion properties of partons, the data sample of central events (0-30% centrality), where a large 
average dijet imbalance is observed, is divided into classes of dijet imbalance. Four Aj selec- 
tions are chosen, which split the central PbPb data sample into approximately equal number of 
dijets: < Aj < 0.13, 0.13 < Aj < 0.24, 0.24 < Aj < 0.35, and Aj > 0.35. For each of these event 
classes, the fragmentation functions are constructed separately for the leading and subleading 
jets. In Fig. |4]ja-d) the fragmentation functions are shown in bins of increasing dijet imbalance, 
from left to right. The corresponding jet pj distributions used in the fragmentation functions 
are shown in Fig Qi-l), illustrating the kinematic range of the measurement and the energy 
imbalance between leading and subleading jets, as selected by the Aj interval. The overlaid 
histogram in the same set of figures shows the pp-based reference distributions. Ratios of PbPb 
data to the pp-based reference are shown in Fig. Qe-h). In the ^ range of 0.5^.0 the PbPb and 
the pp distributions typically agree to within (10-20)% which is smaller than the systematic un- 
certainty in the PbPb measurement, as indicated by the size of the shaded area and open boxes 
for the leading and subleading jet, respectively. Within uncertainties, the PbPb data and pp- 
based reference show the same fragmentation properties, for different jet imbalance in leading, 
as well as subleading, jets. 

6 Conclusions 

The CMS detector has been used to study jet fragmentation properties in pp and PbPb collisions 
at ^ys^ = 2.76 TeV in data samples corresponding to integrated luminosities of about 231 nb^^ 
and 6.8 jih^^, respectively. Jets were reconstructed based on particle-flow objects using the 
anti-A;7- sequential clustering algorithm, with a radius parameter of 0.3. The reconstructed jet 
momenta were corrected to final-state particle level. Dijets were selected consisting of a leading 
jet with px,i > 100 GeV/c and a subleading jet of Pt,2 > 40 GeV/c, with axes that lie within \rj\ < 
2. The azimuthal opening angle A^i 2 between the leading and subleading jet was required 
to be larger than In/ 3. The selected jets were used to construct the high-px component of 
the fragmentation functions by correlating their momentum with the momenta of tracks of 
px > 4GeV/c within a cone of AR < 0.3 around the jet axis. The PbPb results were compared 
to those in a pp-based reference taking into account the different jet momentum distribution 
and the effect of fluctuations in the underlying PbPb event on the jet momentum reconstruction. 
The jet properties were studied as a function of the collision centrality and the dijet transverse 
momentum imbalance. Central PbPb events show a significant excess of unbalanced jet pairs. 
Nevertheless, the fragmentation functions reconstructed in PbPb collisions for different event 
centrality and dijet-px imbalance agree within the measurement uncertainty with the pp-based 
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Figure 4: (a-d) Fragmentation functions for the leading (open circles) and subleading (solid 
points) jets in four regions of Aj in central PbPb collisions compared to the pp reference, (e-h) 
Ratio of each fragmentation function to its pp-based reference. Error bars shown are statistical. 
The systematic uncertainty is represented by hollow boxes (leading jet) or gray boxes (sub- 
leading jet), (i-1) Jet Pt distributions in PbPb collisions in four regions of Aj (not corrected for 
efficiency and not unfolded for pj resolution) compared to the pp-based reference (see text). 
Only statistical uncertainties are shown. 

reference for jets of the same reconstructed momentum. This shows that, after traversing the 
dense strongly interacting medium, partons produced in PbPb collisions are reconstructed as 
jets with a significantly reduced momentum. However, the partition of the smaller momentum 
that remains within the jet cone into high-px particles corresponds to that observed for jets 
fragmenting in vacuum, as seen in pp collisions. 
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